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ABSTRACT 

We study the electron-proton temperature equilibration behind several shocks of the RCW 86 su- 
pernova remnant. To measure the proton temperature, we use published and new optical spectra, all 
from different locations on the remnant. For each location, we determine the electron temperature 
from X-ray spectra, and correct for temperature equilibration between the shock front and the loca- 
tion of the X-ray spectrum. We confirm the result of previous studies that the electron and proton 
temperatures behind shock fronts are consistent with equilibration for slow shocks and deviate for 
faster shocks. However, we can not confirm the previously reported trend of J^/Tp cx l/wg . 

Subject headings: ISM: individual objects (RCW 86 / MSH 14-65) — ISM: supernova remnants — 
radiation mechanisms: thermal — shock waves 



1. INTRODUCTION 

Shocks are ubiquitous in the Universe and occur at var- 
ious scales: from planetary shocks in our Solar System 
to large scale shocks in clusters of galaxies. Astrophysi- 
cal shocks are rather different from shocks on Earth, as 
typical mean-free-paths for particle-particle interactions 
in the interstellar medium are of the order of parsecs, 
which is larger than the scale of the shock transitions 
of interstellar shocks. The fact that we observe sharp 
transitions, implies that the existence of the shock must 
be communicated on scales shorter than this mean-free- 
path. From conservation of mass, momentum and en- 
ergy over the shock front, and assuming an infinite Mach 
number and no cosmic-ray acceleration, 

m = ^-.-^ (1) 

in which Vg denotes the shock velocity , and 7 the equa- 
tion of state of the plasma (7 = 5/3 for a non-relativistic 
gas) . This equation holds for a plasma maximally out of 
thermal equilibrium, implying an electron temperature 
over proton temperature (To/Tp) of the ratio of the elec- 
tron and proton mass (me/mp = 1/1836). For a plasma 
in thermal equilibrium, equation [T] reads: 

2(7-1) 



kT =< kTi >-- 



(2) 



(7+I) 

in which fi is the mean particle mass (0.6 for a fully ion- 
ized plasma with solar abundances). However, if a shock 
is accelerating cosmic rays, a cosmic-ray precursor will 
develop, which heats the incoming medium and lowers 
the effective Mach number of the shock, in which case 



equation [1] or [glare no longer valid 
2009l|Viiik etaL||2010 |. 



see also, Drury et al. 



In most empirical studies on temperature equilibration 
behind supernova remnant shocks, Tq and Tp are deter- 
mined from Balmcr dominated shocks. These shocks are 
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characterized by hydrogen line emission in the B aimer 
series, leading to a line profile consisting of two superim- 
posed lines. These shocks are non-radiative; their energy 
losses in the form of radiation are dynamically negligi- 
ble (this happens for sh ock velocities > 200 kms~^; e.g. 



Draine fc McKee 19931. This implies that the hydro- 
gen lines of these shocks are caused by impact excita- 
tion, not by recombination. The narrow component of 
the spectrum is emitted by neutral hydrogen atoms, ex- 
cited shortly after entering the shock front; its width 
refiects the temperature of the ambient medium. The 
broad component is emitted after charge exchange be- 
tween incoming neutral hydrogen and hot protons behind 
the shock front. The width refiects the temperature of 
the protons behind the shock front. Combining the width 
with the ratio of the flux in the narrow to the flux in the 
broad component, one can determine the shock velocity 
an d hence the e lectron temperature (for a recent review, 
see Heng 2010). This procedure works for shocks which 
have negligible cosmic-ray acceleration. If a shock is an 
efficient accelerator, the cosmic-ray precursor ionizes and 
heats the incoming medium an d therewith increases the 
flux of the narrow component ( Raymond et al.||2011 |. 

The electron temperature for most young supernova 
remnants can also be measured from the X-ray spec- 
trum, as the thermal component of the X-ray spectrum 
is shaped by the electron temperature. Electron tem- 
peratures measured from X-ray spectra are generally not 
the electron temperatures right behind the shock front, 
but rather somewhat further to the inside of the shock, 
where temperatures have already (partially) equilibrated 
through Coulomb interactions. The amount of equili- 
bration is a function of rict: the ionization age, that 
can be determined from the X-ray spectrum simultane- 
ously with the electron temperature, from the ratio of the 
fluxes in different spectral lines from the same species. 



Rakowski et al. (2003) studied shocks of the super- 
nova remnant DEM L71 at both optical and X-ray wave- 
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lengths. The remnant has shock velocities between 555 
and 980 km s~^. Correcting; for te mperature equilibra- 
tion effects, Rakowski et al. (2003) showed that Tc/Tp 
decreases with increasing shock velocity. 

From a theoretical point of view, it is not well de- 
termined how much electron-proton temperature equi- 
libration would to expect behind a shock front. Empir- 
ical studies of shocks with negligible cosmic-ray acceler- 
ation, show that electron and proton temperatures are 
equilibrated for slow shocks {vg < 400 kms~^), whereas 
the degree of t hermal equilibration decreases for faster 
shocks (Rakowski et al. [20031 iGhavamian et al.||2007 

van Adelsberg et al.|2008 
^ 



( Helder et al.||20()9 ). However, its proton temperature is 



71^ 



More specihcally, a relation o: 

has beensuggested for shocks with Vg > 400 



kms""'^ (Ghavamian et al. 20071. This directly implies 
that Te is nearly constant closely behind the shock front 
(i.e. apar t from temperature equih bration effects, equa- 



tion 2 in Ghavamian et al. 



shock velocit y (equation T|). I'he origin of this relation is 



2002), independent of the 



attributed by Ghavamian et al. ( 2007 ) to so-called lower 
hybrid waves ahead of the shock, caus ed by a (mod- 
erate) cosmi c-ray precursor. Note that van Adelsberg 



et al. (2008) could not reproduce Tg/Tp oc v~", using 
pubhshed data of Ha-spectra as well, but with different 
models for interpreting Ha spectra. [Bykov &: Uvarov 
(|l999j) and |Bykov| (|2004| ) point out that for Mach num- 
bers M < ^7m^/m^^ 43 the electrons ahead of the 
shock have a typical thermal velocity exceeding the shock 
velocity, which may lead, through non-resonant interac- 
tions with large-amplitude turbulent fluctuations in the 
shock transition region, to coUisionless electron-heating 
and acceleration. Interestingly, assuming a typical sound 
speed of the ambient medium of 10 km/s, M = 43 corre- 
sponds to a shock speed of ^ 400 km/s. For a cosmic-ray 
accelerating shock it is even more complex, as the incom- 
ing material is already preheated by a cosmic-ray precur- 
sor, leading to a lower Mach number at the main shock. 
This could mean that for cosmic-ray dominated shocks 
the electron-proton temperatures could be equilibrated 
again. 

Here, we investigate the electron-proton temperature 
equilibration behind several shock fronts of RCW 86, us- 
ing X-ray (XMM-Newton/MOS) and optical data. The 
shocks of the supernova remnant RCW 86 are particu- 
larly well suited for this study, as the ionization age of the 
plasm a in parts of the r emnant is low 
(e ' ' 



Vink et al. 



1997 ), this value is low compared to 

the ionization age of most other remnants ^. The low 
value for rict makes the measured electron temperature 
a more direct indication for the electron temperature at 
the shock front. 

Additionally, the shocks of RCW 86 have a large va- 
riety in speeds, as the pr ogenitor exploded in its own 



wind-blown bubble (e.g. Vink et al.||1997 ) and parts 
of the remnant have already hit the dense rim of this 
bubble. As a result, the southwest (SW) rim has slowed 
down; its X-ray emission is mostly thermal and t he shock 
velocity is 500 km s"^ ( [Ghavamian et al.|200T| ). The X- 
ray emission of the northeast (J NE) shock is aomi nated 
by X-ray synchrotron emission (Vink et al. 2006!) , and 
its shock velocity is measured to be 6000 ± 28D0 kms~^ 



low, which is an indication for efficient cosmic-ray accel- 
eration. 

Furthermore, the mostly sub-solar abundances of the 
thermal X-ray emission at the shock front of RCW 86 
indicate emission from shocked ambient medium (rather 



than shocked ejecta, e.g.. 


Vink et al.[ 


19971 [Borkowskil 


et al.[[2001 


Rho et al.[[2002 


1, tying the fitted parameters 
le physics at the shock front. 


of the X-ray spectrum to tl 



2. DATA AND RESULTS 
2.1. VLT/F0RS2 

In this study, we use both published and new opti- 
cal observations. The observation taken from the SW 
of the remn ant was previously published in |Ghavamian| 
et al. (2001) and the obser v ation from the NE rim was 
published in Helder et al. (2009). The parameters of 
the north (N) and east ( E) observations are taken from 
Ghavamian et al. (2007) and the northwest (NW) lo- 
cation from GTiavaniian| ( 1999 ). The exact locations of 
the N and E spectra have been kindly provided to us 
by Parviz Ghavamian (private communica tion). We also 
use pa rameters of a spectrum published in [Long &: Blair 
([1990) , this spectrum is taken from the northern region 
as well. Table 2 summarizes the line widths used in the 
present paper. 

The optical image of Fig. [l] has been obtained with 
VLT/F0RS2 as pre-image for the spectral observation, 
both in 2007 (programme ID 079.D-0735). The image 
is a combination of three exposures of 200s, through a 
narrowband Ha filter (H_Alpha+83). To correct for stel- 
lar light, three similar images through a narrowband Ha 
fiher, shifted 4500 kms"^ (H_Alpha/4500+61) have been 
taken. 

Figure [T] shows the location of the slit, used to obtain 
the spectra. The spectra of the southeast (SE) part of 
the remnant were taken in a single 2733 s exposure, using 
a 1200R grism and a 2.5" slit. This setup resulted in 
a spectral resolution of 350 kms~^, which is sufficient 
for resolving the broad component, however, the narrow 
component of the line probably remains unresolved. The 
data were re duced with standard data reduction steps, as 
described in Helder et al. (2009), resulting in the spectra 
shown in Fig. [5| 

The slit crossed several filaments (Fig. [l} , and we here 
present spectra of the outer and inner nlaments. We 
focus on the line width of the broad component, as this 
is directly related to the post-shock proton temperature. 
We fit the spectra with two Gaussian lines, convolved 
with the resolution of the instrument set-up. Table |2.1| 
lists the line widths obtained from the eastern (SEout) 
and western (SEjn) spectrum in the slit (Fig. [l]). As the 
signal-to-noise of the spectrum of the innermost filament 
is rather low, we used both the inner and middle filament 
for the SEin spectrum. A fit to the innermost filament 
solely revealed a width of the broad line similar to the 
width found for the combined spectrum. 



^ A plasma can be regarded as being in coUisional equilibration 
for net of > 10^^ cm^'^ s. 
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Location 


Broad Ha width 


proton T 


electron T" 






[keV] 


[keV] 


NE 


1100± 60'' 


2.3 ±0.3" 


>1.5 


SEin 


920 ± SO'* 


1.6 ±0.2'' 


0.75 ±0.15 


SEout 


1120 ±40-^ 


2.4 ±0.2'' 


< 0.05 


SW 


562 ± 18"= 


0.60 ±0.04'* 


< 0.05 


N 


325±10 " 


0.2±0.02 


1.95±0.23 


N 


680±70^ 


0.9 ±0.2'* 


1.0±0.2 


E 


640± 35 


0.77± 0.04 


0.78±0.13 


NW 


580± 189 


0.64±0.04 


0.27±0.08 



" At the shock front, derived in this study 
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120091 



>^navamian et~5Z] i[2001| 



Derived in this st udy 
Gha vamian ct al.' |2007[ l 
Cong Ik Blair, |,199r ^ 
i^navamian| | |i9b)t)| l 



2.2. XMM-Newton 

We determined electron temperatures from X-ray spec- 
tra, ta ken with the XMM -Newton/EPIC MOS instru- 
ments ( [Turner et al.||2001|). We chose for the MOS in- 
struments as they have a higher spe ctral resolution than 
the XMM-Newton/EPIC pn CCDs ( [Struder et al.|2001[ ). 
The observations for the E, SEjn and SEout spectrum 
were taken on August 13 (Obsid 0504810201, 75 ks), 
the NE on July 28 (ObsId 0504810101, 117 ks), the N 
and NW on August 25 2007 (ObsId 0504810301, 74 ks 
) and for the SW spectrum on August 23 2007 (.ObsId 
0504810401, 73 ks). The spectra (shown in Fig. |4| were 



extracted at locations close to the regions for which we 
have Ha spectr indicated in Fig. [3] using the XMM 
SAS data reduction software, version 1.52.8. Unfortu- 
nately, the position angle of the XMM telescope was cho- 
sen such that only the EPIC M0S2 instrument observed 
the region of interest for the NE, NW and SW rims. For 
the other spectra, we used both the MOSl and M0S2 
instruments. 

We fitted the extracted spectra with a non-equi librium 
ionization (NEI) model ( Kaastrafe Jansen 1993), com- 
bined with an absorption model, using the SPEX spec- 
tral fitting software version 10.0.0 ( [Kaastra et aI][T996p , 
using the maximum likelihood sta tistic tor Poisson distri- 
butions (C-statistic, |Cash ||1979|). This statistic is more 



appropriate than the classical y'^-method for fitting spec- 
tra w hich contain bins with few counts fWheaton et al. 
19951. For more counts, this statistic asymptotically 
approaches the x^-statistic. For most of the spectra, 
this method was relatively straightforward. Some re- 
gions needed some special attention, as described below. 
The NE spectrum needed an additional power-law com- 
ponent, to account for the synchro tron emission present 
in the spectrum ( jVink et al. |2006 1 . Including the power 
law to account for the synchrotron emission results in 
an unconstrained electron temperature, with a nominal 
electron temperature of 37 keV. However, the continuum 
emission is dominated by synchrotron radiation and the 
temperature is ther efore mostly determined by the very 



weak line emission. 



( 2006] encountered simi- 



Vink et al. 

lar problems but concentrated on a slightly different ex- 
traction region, which had more thermal emission. They 
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Fig. 2. — Ha line in the spectrum of the SEout region (left) and for the SEout region (right). Overplotted in both figures are the the 
individual broad and narrow components (dotted lines) and the total fit (solid line). 





Fig. 3. — RGB image obtained with XMM-Newton. Shown are 
the regions from which the Ha lines have been obtained (white 
circles). The green boxes show regions from which we extracted 
the corresponding X-ray spectra. In this image, red indicates 0.5 
to 1.0 keV, green 1.0-2.0 keV and blue 2.0-4.0 keV. 



reported = 6.7±2.6 keV. The high, but unconstrained 
electron temperature w e found for this st udy, and the 
high value reported by IVink et all (|2006|) suggest that 



the electron and proton temperature are likely close to 
equilibration. However, in or der to be co nservative, we 
take now the 2a lower limit of Vink et aL|2006, (1.5 keV) 
as a lower limit on the electron temperature. 

The SW spectrum can not be fitted adequately with 
a single NEI component (Cash statistic/d.o.f. = 6.4). 
Subsequent analysis with two NEI components with ap- 
proximately solar abundances , which we coupled for the 
components, resulted in a Cash statistic/d.o.f. = 4.5. 



We carried out a different approach carried out by con- 
structing a second NEI component, consisting of pure 
metals. We constructed this component by increasing 
the abundances of O, Ne, Mg, Si, S and Fe with a factor 
of 10^, mimicking a pure metal plasma (a similar proce- 
dure was adopted for the 0519-69.0 remnant in the LMC, 



Kosenko et al. 20101. We found Cash statistic/d.o.f. = 
1.6 for this approach. We found a high Te 



keV) and (1.4 ± 0.1 x 10^" cm' 



(2.2 ±0.3 
■^s) for the pure metal 
plasma component, and a low To and rig for the other 
component (Table [T]). This, together with the complex, 
layered structure oi the remnant in the SW, indicates 
that the pure metal component possibly represents a 
shocked ejecta component in the spectrum and the low 
metallicity component is probably tied to the shocked 
ambient medium. We note that a strong ejecta com- 
ponent with alpha elements was not reported before for 
this remnant, but several studies suggested the presence 
of Fe-rich ejecta giving rise to low ionization Fe-K lines 



famba et aL 2000 Borkowski et al. 2001 Yamaguchi 
al. 2008'). A toliow-up studies with a deeper obser- 



vation IS necessary to shed further light on the emission 
from this region and investigate whether there is indeed 
an ejecta component present in this region. 

The SEout spectrum is fitted with the Fe abundance 
tied to unity, as otherwise this abundance would rise to 
unrealistically high values. We note that releasing this 
component during the fitting procedure will result in an 
even lower electron temperature (0.3 keV). 

Table [l] lists the parameters of the best-fitting spectral 
model. 

3. DISCUSSION 
3.1. XMM-spectra 
The parameters in Table [T] are in general consistent 



with parameters obtain ed in previous work (Rho et al. 
2002|]Vink et al.| |2006| ). In the NE, the dorninant con- 



tribution of synchrotron radiation makes it difficult to 
determine the electron temperature adequately. How- 
ever, we note that the high contribution of non-thermal 
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Fig. 4. — Spectra, corresponding to the regions from Fig. m Red 
lines are the best fitting models. The NE spectrum is muTTiplied 
with 1000, E with 20, SEi^ with 1, SEout with 0.3, N with 0.06, 
SW with 0.0003 and NW with 0.0007. 

X-ray emission to the spectrum makes it difficult to ac- 
curately determine the parameters of the thermal com- 
ponent. Our plasma parameters f or the SW spectrum 
differ from previously determined (Rho et al.||2002), this 



is probably because we fit the spectrum with two com- 
ponents, of which one probably represents the ejecta and 
one the shocked ambient medium, whereas Rho et al. 
used a single component. We choose our region right be- 
hind a Balmer dominated filament, which is non radiative 
and hence has a relatively low density. Additionally, the 
plasma from the ambient medium was probably only re- 
cently shocked. Note that this solution is not necessarily 
unique, but given the high C-statistic values for the al- 
ternative models, we regard this model as an adequate 
representation of the spectrum. The low value for the 
riet for the ambient medium component therefore seems 
appropria te. Alternatively, if we choose our temperature 



( 2002 1 , To and Tp are close to equi- 



similar to Rho et al. 
libration. 

Problems with fitting t he spectrum around 1.2 keV 
were also encountered by Kosenko et al. (20081 for the 
0509-67.5 supernova remnant, 'i'hese deviations are pos- 
sibly caused by uncertainties in the atomic data base 
of SPEX, presumably by the Fe-L line complex. The 
sub-solar abundances in all fits show that we indeed fit 
X-ray spectra from shocked ambient medium as opposed 
to metal-enhanced shocked ejecta. This confirms that 
the measured electron temperatures are related to the 
proton temperatures at the shock front. 

3.2. Proton temperatures 

The width of the broad component is primarily a func- 
tion of the post-shock proton temperature, but modified 
by the velocity dependent cross sections for charge and 
impact excitation. This causes the shape of the broad 
component to deviate from a perfect Gaussian. Nev- 
ertheless this remains a decent approximation (see 
a nd ref erences therein) 



Adelsberg et al. 2008 



oTvan Adelsberg et al 



Fig, 



van 

( |2008[ ) shows that up to 2000 
km/s, the i<'WHM ot the broad line increases linearly 
with the shock velocity. This shows that interpreting 
the line width in terms of proton temperature through 



kT„ 



FWHM/V8 1n2 in kms 
Heni||2010l Table [2l|). 



is a decent approximation (with a 
1 



Rybicki & Lightman 



1979 




10^ 
riet [cm'^s] 

Fig. 5. — Figure with temperature histories for the SEjn location 
(Tp = 1.6 keV), as function rict for different To-values at the shock 
front. From top to bottom, the lines indicate a To at the shock 
front of 1.10, 0.95, 0.80, 0.60, 0.40, 0.20 and 0.05 keV. 



Table |2.1| shows that Tp for SEout is higher than Tp 



for the sEin spectrum, implying a higher shock veloc- 
ity (equation [ij . Note that SEout lies outward of SEin, 
which is suggestive for a higher shock velocity as well. 
For the Northern region of ROW 86, we have the param- 



GhavamianJ 



eters of two Ha spectra ( Long fc Blair|1990 
et al. 2007) at nearly the same location, which appear 
to have significantly different proton temperatures (0.2 
and 0.9 keV respectively). As we do not know which 
value best represents this region, we plotted both values 
for Tp with one corresponding Te in Fig. [6] We note 
that the measured widths of the broad Ha-Imes tend to 
vary a lot (from 32 5 to 543 km/s) in the Northern region 
(Ghavamian 1999), making it unclear to which Tp the 
electron teniperature relates (probably a combination). 
Also, we note that Tp appears t o vary sign ificantly along 
the eastern rim of the remnant ( Ghavamian^,1999) . 



3.3. Temperature equilibration at the shock front 

The proton temperature is obtained from optical spec- 
tra directly behind the shock fronts, as the neutral hy- 
drogen will ionize quickly after entering the shock front. 
To obtain an X-ray spectrum with sufficient signal-over- 
noise, we extracted over a larger region of the remnant 
than the region from which the proton temperature was 
determined. Also, the spatial resolution of the EPIC- 
MOS instruments (6") prevents us from extracting an 
X-ray spectrum from very close to the shock front. 

This means that the electron and proton temperatures 
will have undergone some equilibration at the extraction 
region of the X-ray spectrum, implying that the obtained 
electron temperatures are upper limits (black data points 
in Fig. pi). To compare this to the proton temperatures 
behind the shock front, we need to determine the amount 
of heating the electrons experienced. We assume that 
the elec t rons a re solely heated by Coulomb interactions. 
Spitzer ( 1965 ) derived the coupled differential equation 
for how fast species i and j with different temperatures 
(Ti and Tj) equilibrate as function of time (t): 



dTi _ Tj - Ti 



dt 



tcq (i,j) - 5-87— =2 



T T 



njZfZf\n{A) \Ai A 



For species with atomic numbers Ai and Aj, charge Zj 
and Zj, number density nj and masses mi and toj. A is 
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jiven by 



A 



Solving the coupled differential equation given by equa- 
tion[3]for electrons, protons, He, O, Si and Fe (assuming 
solar abundances), we calculate temperature histories up 
to the measured rict for the measured proton temper- 
ature, for different electron temperatures at the shock 
front (Fig. [s]). This enables us to determine the electron 



temperature at the shock front (Table 2.1 and red data 
points in Fig. pi). 

Equation [l] shows that the proton temperature is re- 
lated to the shock velocity. However, this is based on 
standard assumptions, assuming only a gaseous compo- 
nent. In contrast, supernova remnants have been sug- 
gested to efficiently accelerate p articles which wil l alter 

The 



the standard Hugoniot relations (Vink et al 



presence of X-ray synchrotron and 'i'eV radiation indi 



2010) 



(Bamba et al. 2000| 


Borkowski et al. 


20011 


Vink et al. 


2006 Aharonian et a 


.|2009p, whereas 


the Ha line width 



acceleration must be efficient, wi th a post-shock pre ssure 
contribution of more than 50% ( |Helder et al.]|2009| . 



When a shock accelerates cosmic rays, the cosmic rays 
create a cosmic-ray precursor ahead of the shock, push- 
ing out and pre-heating the incoming medium. This ef- 
fectively lowers the velocity of the material entering the 
main shock, therewith lowering the p ost-shock te mpera- 



ture (Drury et al. 2009 Hughes et al. 2000; Vink et al. 



2010 1. Hence, Vs will not necessarily be the only variable 



to characterize Tp. As we focus on T^-T-p equilibration, 
we plot Te as function of the measured Tp. This is more 
reliable as Vg is a quantity derived through Equation [T] 
which may not be valid in the presence of cosmic-ray 
acceleration ( Helder et al.||2009 ). 

Fig. |6] shows that 7 e ~ 'J p for the data point with 
the lowest proton temperatures. Surprisingly, the North 
Tfi se ems to be very high for the Tp base d on the data 



from Ghavamian et al. (2007). However, Long & Blair 
( 1990p measured approximately the same position and 
found a broader broad component (overplotted in green 
in Fig. [6]). The rightmost data points in Fig. [6] show 
To < Tp, consistent with previous results, however, we 
do not find a relation of To/Tp cx l/wg, which would 
result in a horizontal line, showed as a dashed line (Fig. 

§■ 

4. CONCLUSIONS 

We investigated the proton-electron temperature equi- 
libration behind several shock fronts of RCW 86. The low 
ionization parameter rict of this remnant ensures that the 
current Tc of the plasma is close that at the shock front. 
The electron temperatures were determined from X-ray 
spectra, obtained from XMM-Newton data and corrected 
for equilibration effects, resulting in Fig. [6j We used 
both new and published data to determine proton tem- 



peratures. From this study, we can draw the following 
conclusions: 

- The FWHM of the broad Ha line in the SE of 
RCW 86 are 1120 ± 40 and 920 ± 50 kms^^ for 
the SEout and SEjn spectrum respectively. 
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Fig. 6. — Proton a nd e lectron temperatures {x and j/-axis re- 
spectively) from Table [2?T] and [T] The black data points show the 
electron temperature as calculated from the X-ray spectra. The 
red data points indicate electron temperatures at the shock front 
after correcting for Coulomb equilibration. The dotted line indi- 
cates Tp = Tc. The dashed line indicates the electron temperature 
as calculated for a 400 kms~^ shock with thermal equilibrium be- 
hind the shock front; the relation found by |Ghavamian et al.| | [2007|l 
would be on this line. The points labeled 'N GUY' and 'N LByu 
refer to the same regio n in the North, with Tn based o n parameters 
of , Ghavamian et al._(2007^ and |Long fc Blair] ( [1990^ respectively. 



Tc/Tp ^ 1 for the slow shocks in the E and for 

The N 



the N param eters of Long fc B lair (1990 
parameters of Ghavamian et al. ( |2007) and the SW 
form an exception to this. 

To/Tp < 1 for faster shocks. However, we do 
not find a constant electron temperature for faster 
shocks , as implied by the model of Ghavamian et al. 
( |2007p . Additionally, our results show T^/Tp > 
rrQlrip. 



— To/Tp for the NE region, in which cosmic-ray ac- 
celeration appears to be efficient, is not well deter- 
mined. However, we have found evidence for rel- 
atively high values of To, suggesting that Te/Tp is 
close to one is this region. This may be attributed 
to efficient cosmic-ray acceleration, as this tends to 
lower the Mach number of the main shock. 
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